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ABSTRACT 

When tracking a target particle that is interacting with nearest neighbors in a known way, positional data of 
the neighbors can be used to improve the state estimate. Effects of the accuracy of such positional data on the 
target track accuracy are investigated in this paper, in the context of dusty plasmas. In kinematic simulations, 
notable improvement in the target track accuracy was found when including all nearest neighbors in the state 
estimation filter and tracking algorithm, whereas the track accuracy was not significantly improved by higher- 
accuracy measurement techniques. The state estimation algorithm, involving an extended Kalman filter, was 
shown to either remove or significantly reduce errors due to "pixel- locking" . It is concluded that the significant 
extra complexity and computational expense to achieve these relatively small improvements are likely to be 
unwarranted for many situations. For the purposes of determining the precise particle locations, it is concluded 
that the simplified state estimation algorithm can be a viable alternative to using more computationally-intensive 
measurement techniques. 
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1. INTRODUCTION 

Dusty or complex, plasmas consist of a low-density neutral/ion/electron plasma containing suspended "dust" - 
negatively-charged macroscopic particles. Dusty plasmas occur in space as well as in various terrestrial discharges 
ranging from lightning to industrial applications. 1 In more controlled environments, dusty plasmas offer a unique 
testbed for exploring very complex and fascinating physical processes on a kinematic level. In all cases, precise 
determination of properties of the dusty plasma is a top priority. In particular, the industrial applications of 
dusty plasmas could benefit from being able to determine, and subsequently control, physical properties of a 
dusty plasma. 

Many properties of a dusty plasma can be inferred from the behavior of the dust particles, which can be 
non-invasively observed over time via laser illumination and a digital camera, as in Figure [I] which is used 
in experiments on Mach cones 2 and shock- waves. 3 The particles acquire a charge due to collisions with ions 
and electrons in the plasma, causing them to interact with each other and the plasma, with often fascinating 
results. For example, the dust can form crystalline structures or behave in a liquid-like phase (or exhibit phase 
transitions^), on timescales amenable to observation at the individual particle level. The precision of one's 
knowledge of the dust behavior is crucial to the reliability of any conclusions made about physical properties of 
the system)^ There are two ways to improve the precision of determining the dynamics of individual particles: 
more sophisticated measurements, and information processing. The performance of each are considered in this 
work. 

Determining the tracks of multiple, interacting targets is a non-trivial task. Dust particles are typically 
observed using a video camera, their locations revealed using computer-based image-processing techniques, and 
their tracks (connected positions in time) formed using track-association algorithms. Such techniques and algo- 
rithms are available in varying levels of sophistication, with commensurate levels of accuracy, precision, resource 
consumption, etc. That is, more sophisticated techniques typically provide better precision and consume more 
resources such as time and computer memory. This article considers tracking performance and, specifically, a 
trade-off between resources and accuracy/precision. Motivation for such considerations comes from the ulti- 
mate goal of implementing closed- loop control for such experiments, for which the tracking procedure must be 
automated. 



Previous work on tracking individual dusty plasma particles has been somewhat unsophisticated. A common 
technique is particle tracking velocimetry (PTV), where the distance and direction of particle travel is obtained 
from consecutive images and used to aid in linking particle locations into tracks. 5 6 This process can be auto- 
mated, but more can be done. A higher level in sophistication for particle tracking is offered by the Kalman 
filter JZE which combines a prediction model with measurements to obtain a recursive state estimation algorithm 
that can be optimal in the minimum-mean-square-error sense for linear dynamics. A simplified Kalman filter has 
been used on simulated data in the context of a dusty plasma crystal,^ where it was shown to perform very well 
in many situations. Slightly more sophisticated still is the state-estimation algorithm known as the Extended 
Kalman filter (EKF) J3H which is considered in this work and will be introduced in Section K3j 




Figure 1. Schematic of the dusty plasma experimental setupP^ The dusty plasma is contained wholly within a chamber 
as shown in (a), with the dust illuminated by a laser sheet and imaged from above by a digital camera as shown in (b). 



2. SIMULATED DYNAMICS 

Simulating experiments affords this study two primary benefits. It allows for a larger exploration of "parameter 
space" than might be achievable experimentally in a practical period of time. In addition, it allows for the 
accuracy of the algorithms to be quantified because the "true" particle locations are known. In the future these 
algorithms will be developed for and tested on experimental data, so it is important to be as realistic as possible 
at this early stage. To this end, the parameters and images of recent dusty plasma experiments^!!^ are used as a 
guide for this work. 

Three characteristic length scales of importance for dusty plasmas are the particle radius R, the Debye screen- 
ing length Ad and the inter-particle separation r. For the case where R <C Ad < r, the dust particle interactions 
can be described welJ^ by a point-charge approximation. Each particle of charge Qd experiences forces from 
other particles directly, as well as indirectly via the plasma (such as drag forces due to ion motionmotion of 
neutrals within the plasma). Under the point-charge approximation, the effective potential experienced by one 
dust particle (labeled j) due to another (labeled k) consists of a screened (exponentially decreasing), repulsive 
Coulomb interaction of the Debye-Hiickel/ Yukawa formP^l 
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Here fco = V(47T£o) * s Coulomb's constant, the plasma screening distance is the Debye length (Ad), and the 
particles are separated by f 3 ^ = rr j,k (note that this separation is time- varying) . Assuming a circular two- 
dimensional geometry (no angular variation), and defining f = r/Ao, yields the effective force between two 
interacting dust particles (the gradient of the potential) to be 

which is zero at a particle separation of r — (1 + y/3) Ad = ro- For r > ro, two simulated dust particles attract 
each other weakly, and for r < ro they repel strongly. When used as the initial separation between simulated 
dust particles, small (thermal, etc.) deviations can be expected about this minimum, resulting in a crystal-like 
structure for which an EKF is used for the state estimation procedure (see Section [3| . When the combined 



effect of surrounding particles and a global confinement potential are included, dusty plasma particles tend to 
align themselves in a hexagonal lattice,^ about which small (thermal, etc.) deviations can be expected after a 
transient period. This situation is considered here. 

The total force on the jth particle is then the sum over k for the two-particle interactions ([2| and the particle 
dynamics are determined using a fourth-order Runge-Kutta numerical integration of the equations of motion for 
position q(t) and velocity v(t): 
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W(t))=0, 
W(t)W(t')) = S(t-t'). 



Brownian motion of the dust particles due to collisions with plasma ions has been included in a standard way: 
a zero- mean delta-correlated stationary Gaussian process W(t) = dW{t)/dt, which satisfies 

(5) 
(6) 

The magnitude of a particle's position fluctuations induced by Brownian motion is related to the local temperature 
T by the Einstein relation: D = k^T '/jm, where /cb is Boltzmann's constant, 7 is a particle damping coefficient, 
and m is the particle mass (~ 5 x 10 -13 kg in Ref. |9j). 

2.1 Images 

In a dusty plasma experiment, observations consist of a time sequence of images taken (typically) at regular 
intervals of At. Thus, when developing a state estimation and tracking algorithm using simulated data, it is 
beneficial to reproduce images that approximate those obtained experimentally. 

Experimental images from Ref s. [3l [l2| and [13] provide the benchmark for the images generated from simulations 
in this work. The raw images from the experiments reported in Ref. 13 compare very well to simulated images, 



as shown in Figures 2(a) and 2(b) 




Figure 2. Zoomed versions of (a) 



(c) (d) 

captured image from experiment, (b) a simulated image without camera defocussing, (c) 
a simulated single hexagonal lattice "cell" consisting of a central target particle surrounded by the six nearest neighbors, 
(d) a simulated image with simulated camera defocussing to increase pixel-resolution and reduce/remove pixel- locking 
errors. Images from experiments and simulations compare very well. 



2.2 Nearest-Neighbor Approximation 

Given the exponential suppression of the Yukawa force in Equation it is reasonable to consider the relative 
effect of the Nearest-Neighbours (NNs) and the Next-Nearest-Neighbors (NNNs) on a "target" particle. For a 
crystal-like state where r > Ad, the exponential suppression factor suggests that the NNNs have relatively little 



influence on a target's dynamics. We will now show that this is true by analyzing propagation of uncertain- 
ties. Uncertainties in the two-particle separation r propagate into F(r) as SF(r) = \{dF/dr) \ Sr, where dF/dr 
calculated from Equation Q is 
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and the uncertainty in the particle separation will typically be Sr < V25q, where Sq is the largest measurement 
uncertainty in either dimension. That is, the average error in locating a particle's center. In this work, Sq ~ 
0.1 pixels. Consider now the ratio between uncertainties in the two-particle forces due to NNs and NNNs: 
£F(tnnn)/^(7nn)- If this ratio is sufficiently small, then we can safely omit the NNNs from the tracking 
algorithm. Defining = tnnn/^nn > 1 and working in units of tnn = ^nn/Ad, we find that 
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which is approaching insignificance (f(0) < 10%, say) for larger than approximately 2, as shown in Figure [3j 
Thus, for crystal-like states where tnnn ~ 2tnn, the tracking algorithm used need only consider the NNs and 
not the NNNs. This is a great reduction in complexity. 




Figure 3. Plot of Equation ([8| showing the effect of considering the next nearest neighbors (NNNs) in the tracking 
algorithm filter. The uncertainty in the predicted force on the tracked central particle is insignificant (f(0) < 10%, say) 
for >~ 2, which is always true for crystal-like structures of dusty plasmas where tnnn ~ 2rNN- If the NNNs move 
much closer than this, say < 1.5, then the influence of the NNNs may need to be considered in the tracking and state 
estimation algorithm. 



2.3 Simulating a single hexagonal lattice cell 

As shown in the preceding section, for a dusty plasma in a crystal-like phase, a target particle is only signifi- 
cantly affected by its NNs. Therefore, for the purposes of ascertaining (and quantifying) how the knowledge of 
surrounding particles can benefit the tracking process, this work focusses on a single cell of a hexagonal lattice, 
as seen in the experimental and theoretical work of Refs. [2| [3J [5| [9| and |T2j for example. A small set of seven 
particles was simulated, representing a single cell of NNs surrounding a target particle in a hexagonal lattice as 
can be seen in Figures [2(c)1 and [2(d)] The seven particles are initially positioned at an equidistant separation 
near potential minima (total potential: Yukawa plus confining potential) and interact via the Yukawa-type of 
potential ([T]) as described in Section [5] The shallow, parabolic, global confining potential is used to mimic the 
combined effect of many other surrounding particles and an experimental confining potential P These assump- 
tions produce a reasonable model for a small subset of a larger two-dimensional dusty plasma in a crystal-like 
state. 



3. STATE ESTIMATION AND TRACKING 



3.1 Measurement: image-processing 

Once images are obtained from a dusty plasma experiment (real or simulated), they are processed to determine 
the centroids (sub-pixel resolution location) of each particle. Herein, these centroids will be referred to as the 
particle measurements. The centroid location is achieved using one of a number of techniques described in 
Ref. [14] The simplest of these is the Threshold Method (TM) where contiguous pixels above a certain threshold 
are considered to be a particle, with the pixel center-of-mass used for the centroid. The Moment Method (MM) 
works in a similar way to the TM, but calculates a center-of-mass weighted by the particle's pixel-intensities. 
A third technique known as the Moment Method with Filter (MMF) involves applying a spatial bandpass filter 
(typically Gaussian) for smoothing and background subtraction prior to applying the MM to find the weighted 
center-of-mass. In typical order of increasing precision, these techniques are ranked TM, MM, MMF. There 
exist more sophisticated (and often more precise) measurement techniques, such as the Linear Quadratic Kernel 
Method™ for which the increased computational expense cannot necessarily be justified for all applications, 
particularly if combining the measurements with a state-estimation algorithm that can improve the precision of 
the estimated positions for less computational expense. 

Another option for improving the precision of dusty plasma particle measurements involves a hardware 
modification: defocussing the camera lens to obtain a "smeared" image to which the simpler (and faster) TM and 
MM methods can be applied. This approach (used here) helps to avoid errors due to "pixel-locking"E3El! where 
certain sub-pixel locations are preferred due to low pixel-resolution caused when a particle only illuminates a small 
number of pixels. The particles in Figure 2(c) are of the order of 3 pixels across, whereas those in Figure 2(d) are 
approximately twice as large. The level of smearing can be modified to obtain a desired measurement precision, 
but it should be sufficiently moderate so as to not affect the measurement-track association process discussed in 
Section IO 



3.2 Measurement-track association 

Each measurement must be either associated to a track, dismissed as a false alarm, or saved in memory as a 
possible new track. Along with the task of allowing for missed detections, this decision process is known as 
measurement-track association. In multi-target situations this can be a highly nontrivial task, especially when 
two or more targets are close to one another, for example. For simple situations where the targets remain spatially 
well-separated between subsequent measurements (Euclidean separation is greater than the measurement and 
prediction noise) , a simple radar-like technique can be used to locate the measurement nearest to each track. This 
is the case for crystal-like dusty plasma structures close to equilibrium. This rather simple scenario provides a 
benchmark for measuring the performance of more sophisticated measurement-track association algorithms that 
are necessary for highly nonlinear dynamics such as Mach cones^ and shock- waves P 

3.3 The Extended Kalman Filter 

The state estimation for the particle tracks uses a standard discrete-time Extended Kalman Filter (EKF)EE 
The state of each particle being tracked is represented by six variables - two positions g x and q y , two velocities 
v x and v y , and two accelerations a x and a y - at each time step k = 0,1,2 — The state is contained in a 
six-element vector: x(k) = (<? x , u x , a x , g y , v y , a y ) T and the measurement for each particle is represented by a 2 x 6 
measurement matrix H: 

1 
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so that the expected measurement for each track, given the previous measurements up to and including k — 1, 
is given by 

z(k\k -1) = H- x(k) = ( ) . (10) 



q y (k) 

The expected and actual measurement (z(k)) at time step k are combined to form the innovation (the difference 
between the expected positions and the measured positions) and added to the state vector as in the standard 
Kalman filter, 

x(k+) = x(k) + K(k) - (z(k) - z(k\k)), (11) 



where x(k+) indicates the updated state vector, the Kalman gain is given by K(k) = S x (k) ■ H T ■ (H ■ S x (k) ■ H T + 
-R(fc)) -1 , S x (k) is the estimated covariance matrix for the track states, and R(k) is the expected covariance matrix 
for the measurement process being used. 

After each time step, the particle state vector is predicted forward to the next time step using a discretization 
of the dynamical equations used in the simulation and given in Section [2] 
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where At is the time step. The corresponding equation for the predicted covariance is 

S x (k + 1) = f(k) ■ S x (k + ) ■ f'(kf + Q(k), 



where Q(k) is the process noise covariance, and 



(13) 



(14) 



x=x(kj r ) 



This constitutes a piecewise linearization of the nonlinear particle dynamics represented by the simulation. This 
linear approximation should be accurate as long as the errors in the state estimates are small compared to the 
nonlinear nature of the underlying potential function. For the cases considered here, with the simulated dust 
particles in a crystal-like state, the relative movement of the particles around their equilibrium positions will 
be small enough for this approximation to hold. In more dynamic dusty plasmas, checks would be required to 
ensure that the estimated errors were not exceeding the ranges required for the perturbative expansion of the 
inter-particle forces about the current state estimates. 

3.4 Process noise tuning 

A tunable parameter for any EKF (or simple KF) is the process noise magnitude <tq, which models the standard 
deviation for the highest-order time-derivative of position that is not included in the prediction model - in this 
case, the jerk (the time-derivative of acceleration). As these typical fluctuations are generally unknown a priori, 
<tq becomes a design parameter that needs to be optimized to help minimize modeling errors. This "process 
noise tuning" involves minimizing the average mean-square or root-mean-square (RMS) error as a function of <tq. 
This is one example of where virtual experiments are beneficial because many experiments (perhaps hundreds or 
thousands) can be simulated to obtain good averages. For a typical simulation ensemble of size N = 1000, <jq 
was varied over a finite range and the resulting average RMS error in the state estimate of the target particle's 
position, (^ x ,%) was calculated ad^ 
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where Aq x ^(k) is the deviation of the estimated x-position from the true x-position at time t{k) (and similarly 
for the y-direction) . The long-time average of ( p~5| ) provides a convenient measure of the filter's performance as 
a function of <tq. This is shown in Figure [4] (note the logarithmic horizontal axis), where data corresponding 
to TM measurements is blue and MM measurements is red (and bold). For increasingly large process noise 
magnitude, confidence in the prediction model decreases so that the measurements are "trusted" more by the 
filter and the RMS error in the estimated state asymptotically approaches the measurement noise level. An 
optimal process noise magnitude yields a minimum RMS error in the estimated state - note that the TM state 
estimate for a tuned filter performs better than the (non-optimized) MM measurement level, providing very 



encouraging support for a software-based approach to improving particle-locating precision, as opposed to what 
could be a laborious hardware optimization process in a real experiment in order to obtain comparable errors 
sizes. In any case, the state estimate outperforms the measurements for sufficiently large process noise, showing 
that it's almost always beneficial to use any knowledge you have of the target particle dynamics to improve track 
precision via an EKF (the payoff decreases exponentially for less certain knowledge as the information processing 
costs outweigh the reduced improvement in precision). 

The optimal values of process noise magnitude were crQ 5 TM(At) 2 = 0.05 pixels per frame and crQ 5 MM(At) 2 = 
0.1 pixels per frame, where the sampling period used corresponds to imaging the virtual experiment at 200 
frames per second: At = 0.005s. The lower precision of TM measurements should result in a higher level of 
confidence (smaller ctq) in the prediction model, and so the RMS errors in the state estimate are minimized for 
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Figure 4. For a single ensemble of 1000 virtual (simulated) experiments, multiple EKFs were run using a different 
process noise magnitude, expressed here in units of pixels per frame: ctq x (At) 2 . The optimal values (corresponding 
to fluctuations in particle positions of one-tenth and one- twentieth of a pixel per frame for MM and TM measurements, 
respectively) minimize the average RMS error (given here in pixels) over long times. 




3.5 Filter consistency 

For static parameter estimation, consistency of an estimator is defined as convergence of the estimate x to the 
true value x. For dynamic parameter estimation, estimator consistency is defined in terms of conditions on the 
first two moments of the estimate: 7 

E[x-x) = 0, (16a) 
E[{x = x}{x = x}'\ = P. (16b) 

The first is satisfied by an unbiased estimator (zero-mean estimation error). The second condition is that of 
covariance matching between the MSE of the estimates (left-hand side) and the filter-calculated covariance 



(right-hand side). The two conditions in Equation (16) can be tested simultaneously using a x-squared test, as 



described in many undergraduate-level statistics texts - Ref. u\ describes it in the context of a Kalman filter. 
When the test statistics reside within an acceptance window, then the filter is determined to be consistent, which 
was the case in this work (although they are not presented here). 



4. RESULTS AND DISCUSSION 
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(a) Sample particle trajectory 



(b) Optimal RMS errors over time 

Figure 5. Shown in |(a)| is a sample trajectory showing the true evolution of the target's Y-position (thin black line), a 
spread of TM (blue +) and MM (red x) measurements and the corresponding state estimates for the tuned process noise 
values given in Figure [4] Note the significantly larger deviation of the measurements compared to the estimates. Shown 
in (b) are the RMS errors in measured and estimated positions as a function of time, again for the tuned values of process 



noise. The sample period was At = 0.005s (200 frames per second). 



Figure [5(a)] is a sample trajectory from the early stages of a single virtual experiment. The true Y-position 
of the central particle is shown as a thin black line, with the MM and TM measurements shown as crosses (x) 
and plus signs (+), respectively. Shown are the state estimates generated by the EKF after processing each 
measurement and weighing it with the prediction model - the red line (MM) and blue dashed line (TM). Note 
the transient period before the filter "settles down" and tracks the true state quite well. 

Figure |5(b)1 shows the RMS errors for the measurements (dotted lines) and estimated positions (solid lines) of 
the target particle when tracking only the target, but feeding the NN measurements into the prediction stage of 
the EKF. When tracking all 7 particles in the hexagonal "crystal" cell, the increased precision in the tracked NN 
positions did not significantly improve the target particle track. This is reasonable because the measurements 
were quite good, as shown earlier. For situations where the measurements are not so good (due to noisier images, 
for example), or the NN particles are moving faster, tracking the NNs (as opposed to just measuring them, which 
can represent a massive computational saving when tracking many particles) is expected to produce noticeably 
better results for the target track, albeit at a noticeably higher computational expense due to the increased 
number of states in the filter. For longer times than shown, the RMS error in the estimated position exhibits 
low-magnitude oscillations at a rate corresponding to the rate of oscillation of the interacting dust particles about 
the initial positions. This explains the apparent slight increase in the state estimate RMS errors at later times 
in Figure [5(b)] 



As shown in Section 3.4 , the process noise magnitude ctq can be tuned to improve the EKF performance 
(minimize the errors in the estimated state). In the same way that the magnitude <jq represents confidence in 
the prediction model, the measurement noise magnitude <jr represents confidence in the measurements. Figure [4] 
shows that the actual measurement errors were quite small: <tr 5 tm ~ 0.16 pixels, and ctr 5 mm ~ 0.1 pixels. In 
this work, no knowledge of the measurement precision was assumed and the so-called 4 "worst-case scenario" 
of <jj^ = 0.1 was used. Overestimating the measurement noise in this way is usually not as detrimental to the 
resulting performance of a Kalman filter as underestimating it, as shown for a simplified Kalman filter in Ref. [4j 
Therefore, it is expected that using a measurement noise size in the EKF that is closer to the true value would 
not significantly improve the EKF performance. Note that a similar tuning process can be performed for the 



measurement noise as was done for the process noise in Section 3.4 
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Figure 6. Sub-pixel maps for the central particle's (a) true positions, (b) TM-measured positions, (c) MM-measured 
positions, (d) TM-based estimated positions, (e) MM-based estimated positions, across 1000 virtual experiments. Note 



the presence of pixel- locking in the TM measurements in (b) (and to a lesser extent in the MM measurements in (c) ), and 
its subsequent removal by the EKF in (d) 



There are a few notable potential sources of error in tracking particles in dusty plasma experiments (both 
virtual and real). They can be categorized as measurement-based, model-based and numerical. The latter 
come in the form of round-off errors that can be minimized by using the Joseph form of the covariance update 
equation in the EKFP Modeling errors, where an incorrect model is used for the system dynamics, can lead to 
incorrect predictions, and over- or under-confidence in predictions (versus measurements) due to sub-optimal 
process noise values - overcome here by the process- noise-tuning process. The most serious measurement-based 
errors for the rather unsophisticated techniques considered here (TM and MM) come in the form of "pixel- 
locking" JHE3 resulting from insufficient pixel-resolution for the particles. This results in certain subpixel regions 
being preferred, as shown in Figure [6] where the measured subpixel locations (Figure [6 (b) | and [(c)] ) of the target 
particle in 401 images from each of 1000 virtual experiments (401,000 data points) are plotted alongside the true 
) and estimated (Figure |6(d)| and |(e)[ ) subpixel locations. Pixel-locking was reduced for the MM 
measurements (see Figure |6(c)[ ) in a fairly standard way - "smearing" the particles slightly by defocussing the 
virtual camera. Remarkably, although this technique was completely unsuccessful for the TM measurements (see 
Figure [6(b)] ) , the filter almost completely removed the pixel- locking effect as shown in the subpixel map for the 
TM-based estimates in Figure [6(d)] These maps were representative of all the particles. They show that strong 
pixel- locking occurred for the TM measurements (and to a lesser extent for the MM measurements), but was 
removed after processing these measurements through the filter. 



5. CONCLUSION 



An extended Kalman filter (EKF) algorithm has been designed and implemented for tracking interacting dusty 
plasma particles. The algorithm was tested in virtual experiments (simulations) which allowed for quantification 
of performance aspects such as the precision of locating particle positions. Attention was paid to generating 
images that closely resemble real experiments. The algorithm was tested for a crystal-like phase of the dusty 
plasma. This simple scenario provides a benchmark for measuring the performance of more sophisticated al- 
gorithms that will be necessary for treating more complex dynamics. An analytical result showed the validity 
of an approximation to only nearest-neighbor interactions for this situation. Starting with relatively imprecise 
measurements of particle positions using the threshold and moment methods, it was shown that the filter algo- 
rithm improved the precision, as well as removing or significantly reducing errors due to "pixel- locking" . The 
RMS values of position errors exceeded 0.1 pixels for the measurements, which were improved by the filter to 
approximately 0.06 pixels - approaching the precision obtained in Ref. |l4| (« 0.03 pixels) where a much more 
sophisticated and time-consuming measurement technique was used (local quadratic kernel method). This sug- 
gests that when such "cheap" measurements are used in conjunction with an EKF, they may be sufficient for 
obtaining excellent sub-pixel precision for particle positions - a task that was previously achieved through hard- 
ware modifications and/or the aforementioned sophisticated measurement (image-processing) techniques. This 
particular result bodes well for eventually applying such a particle-tracking filter in real-time, perhaps as part of 
a feedback loop. 

It is expected that the algorithm's will perform quite wellperformance in experimental situations will be 
comparable to the results presented here, but performance will degrade the further a dusty plasma deviates 
from a crystal-like state. Future work includes applying the algorithm to experimental images, extending the 
algorithm to cope with highly nonlinear dynamics such as shock-waves and Mach cones, and considering more 
advanced state estimation and measurement-track association techniques. The present results will provide a 
useful performance benchmark for future work. 
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